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Abstract Limit theorems for the time average of some 
observation functions in an infinite measure dynamical 
system are studied. It is known that intermittent phe- 
nomena, such as the Rayleigh-Bcnard convection and 
Belousov-Zhabotinsky reaction, are described by infinite 
measure dynamical systems. We show that the time av- 
erage of the observation function which is not the L^{m) 
function, whose average with respect to the invariant 
measure m is finite, converges to the generalized arc- 
sine distribution. This result leads to the novel view that 
the correlation function is intrinsically random and does 
not decay. Moreover, it is also numerically shown that 
the time average of the observation function converges 
to the stable distribution when the observation function 
has the infinite mean. 

Keywords Non-stationary Chaos • Infinite measure • 
Generalized Arcsine Law • Non-equilibrium state 



1 Introduction 

Recently, the 1// power spectrum and the power law 
phenomena, which are closely related to the intermittent 
phenomena [T], have been studied in various systems. 
Examples are the 1// power spectrum in the Rayleigh- 
Benard convection [2], Belousov-Zhabotinsky reaction [3j, 
fluorescence intermittency in single nanocrystals [4] and 
the power law decay of the earthquake phenomena [SI 
[6]. It is also known that such power law phenomena 
are clearly observed in Hamiltonian systems and non- 
hyperbolic dynamical systems [HIH]. It is a remarkable 
problem that the theoretical meaning of the 1// power 
spectrum has not been completely elucidated. 
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In ergodic theory the time average can be replaced 
by the ensemble average. However, in the intermittent 
phenomena, this replacement is not always guaranteed. 
Actually, it has been pointed out that the time average 
of some observation shows anomalous behaviour, that is, 
the time average does not converge to a constant value 
and becomes intrinsically random [QlllOj. Therefore it is 
important to analyze the behaviour of the time aver- 
age in the intermittent phenomena from the viewpoint 
of ergodic theory. Infinite measure dynamical systems 
are examples of dynamical systems describing such in- 
termittent phenomena. The recent study of infinite er- 
godic theory tells us that the time average of some ob- 
servation functions converges in distribution [TTlfT^IT^ 
[T4] . For example, the scaled time average of the L^(m) 

function g{x), X]fc=o ffC^'^O/'^nj converges to the Mittag- 
LefSer distribution, where a„ is the proper sequence. 

The purpose of this paper is to study the time aver- 
age of the non-L^(TO) function in infinite measure sys- 
tems. In the previous works, it is already known that 
the distribution of the occupation time of some interval 
in infinite measure systems converges to the generalized 
arcsine distribution |151ll6j . In this paper we clear the 
class of the observation function whose time average con- 
verges to the generalized arcsine distribution using the 
modified Bernoulli map which is the typical example of 
a one-dimensional infinite measure dynamical system on 
[0,1]. 

The paper is organized as follows. In §2 we review 
the modified Bernoulli map. In §3 we show the distribu- 
tion for the time average of the i;\,c(0, 1) function with 
finite mean using a modified Bernoulli map. Theorem 1 
and Theorem 2 are the generalized arcsine laws, which 
are well known in random process, in the dynamical sys- 
tem. Theorem 3, which is the main result in this paper, 
generalizes the observation function. In §4 our results are 
applied to the correlation functions. In §5 we numerically 
demonstrate the distribution of for the time average of 
the ijocC'^' 1) function with infinite mean. §6 is devoted 
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to summary and an approach toward the ergodic prob- 
lems of non-equihbrium statistical mechanics. 



2 Reviews of the modified Bernoulli map 

In this section we review the statistical properties of the 
modified Bernoulli map. The modified Bernoulli map is 
one dimensional map on [0, 1] defined by 




[l-XnY 



Xn e /o = [0, 1/2] 

a-„e/i = (1/2,1]. 



This map has two indifferent fixed points (x = and 1), 
at which the invariant density p{x) is not bounded. The 
invariant density p{x) can be written as 



p{x) 



1~B 



(2) 



for a; and x ~ 1. Therefore for B > 2 the invari- 
ant density cannot be normalized, that is, the invariant 
measure becomes the infinite one [7l [T7lll8| . 

It is important that renewal processes are constructed 
by the sequences of the modified Bernoulli map [Hj . Ac- 
tually, using the symbolic sequence (t„ = cr(a;„), where 
<t{x) = —1 {x G Iq) and (t(x) ~ 1 (x € Ii), one can 
define the renewal when the value of cr„ changes, namely 
fnCTri+i = ^1- As the time interval between r — 1th 
and rth renewal denoted by (r > 2) is an inde- 
pendently identically distributed random variable, and 
this probability density function (p.d.f.) f{n) is given by 
/(n) oc (n - 1)-^ (/? = B/{B - 1)) for n » 1 |20], this 
process is regarded as a renewal process. In what follows, 
we consider the time interval between successive renewals 
as continuous random variables, i.e., p.d.f. is given by 



fix)^{(3~l){x-l)-''' U3 = B/{B-1)). 



(3) 



It is noted that the p.d.f. of the first renewal time Xi de- 
pends on the initial ensemble of the modified Bernoulli 
map pni^ . When the initial ensemble is the invariant 
density for the first passage map T"*^^^(a;) with respect 
to E = [ei, 62], whose endpoints are the solutions of the 
equation Tx ~ 1/2 for ei < 62 and n(x) = 1 + min{n > 
: r"(a;) G E}, the p.d.f. of Xi is same as Eq. ®, 
namely the ordinary renewal process. However, when the 
initial ensemble is the invariant density for the modified 
Bernoulli map, the p.d.f. of the first renewal time Xi 
is given by fi{x) = (1 — F{x))/fi, which is completely 
different from Eq. ([3]), where F{x) is the cumulative dis- 
tribution function of f{x) and fi is the mean value of X^ 
[T9] . In the previous papers pT1[22] we clearly demon- 
strate the dependence of the statistical laws, namely, the 
renewal function and the correlation function, on the ini- 
tial ensemble. 



3 Generalized arcsine law 

Firstly, we analyze the behaviour of the time average of 
the following function: 



i{x) = c - ?! 



(4) 



where a,b G R \ {—00, 00}^ It is noted that I{x) is not 
the L^{m) function. 

We review Lamperti's generalized arcsine law for the 
modified Bernoulli map 

Theorem 1 (Lamperti's generalized arcsine law) 

Let X„ he the time interval between the successive re- 
newals in the renewal process constructed by the modified 
Bernoulli map and Nn be the occupation time in Iq, that 
is, Nn = Xi -I- X3 -I- ... -I- X„ when xq G Iq and n is odd. 
Then 

lim Pr(iV„/n < x) = Ga{x) (5) 

n — >oo 

exists, where a = /? — 1 and the p.d.f. G'^{x) is given by 
sinTra x'^ {I - x)"'^ + x"''^ {I - x)" 



G'^ix) 



IT -I- 2a;"(l - a:)" cosTTo; -I- (1 - 



(6) 



The distribution Ga{x) is called the generalized arc- 
sine distribution. 

Proof. In [24] the limit probability (O exists if and 
only if there exist constants c and a such that 

lim EiNn/n) = c (7) 

n — >oc 

and 

{l-x)u'{x) 

lim — j r-r- = (8) 

where E{.) is the ensemble average with respect to the 
initial density and u{x) = f{'n)x^ is the generating 

function of the f{n). In what follows, we check the above 
conditions. 

First, we define the average of the function l[o,i/2) {T"'x) 

as 

E,, = E{l[o,i/2){T"x)). (9) 

According to the initial density po{x) converges to 
the invariant density under proper normalization: 

1 



WnP''po{x) 



p{x) as ri — > 00, (10) 



^r(a)r(2-a), 
where P is the Perron-Frobenius operator and 



Wn 



logn (S = 2) 
{B > 2), 



,1-a 



(11) 



where a = (3—l. The invariant density p{x) is symmetric 
with respect to the axis x = 1/2. Therefore 



E„ — > — as n — > 00, 
2 



(12) 



^ The special case, namel y a — —1 and 6 = 1, is shown by 
using the renewal theory in |23) . 
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and 

E{N^ln) 



1 " 
n ^-^ 

k=l 



— a.s n 
2 



oo. 



(13) 



By Karamata's Tauberian theorem it is easily con- 
firmed that the condition ([8|) holds when a = P — 1: 



hm 

x^i^ 1-F{x) 



1. 



(14) 



The p.di. oiNr,/n in [24] under c 1/2 and a = /3-1 
implies the p.d.f. □ 



Remark 1 In the case of the modified Bernoulli map 
the measure of the sets [0, e] and [I — e, 1] are not finite. 
Therefore we exclude these sets in Eqs. (| and 
When the measure of the other sets are not finite, Equa- 
tions and (ji30p must be changed to exclude these 
sets. 

Theorem 3 Let g{x) be the Lj^^{0, 1) function with fi- 
nite mean w.r.t. the invariant measure m and ^(0) = 
a,g{\) ~ h. Further, there exists 5 such that < 5 < 1 



. mi m 1 ii 1- i -1 i- r o.'Ti'd qix) IS continuous m 0, i5 U 1 — 15, 1 . Then the time 

Usma; ihcorcmlll one can know the distribution oi i / \ ■ l j_ -l j_ j_ i 

, . \ ■ J- i 1 averaqe ojaix) converqes m distribution to Ya a b, where 

the time average OI 7 (x) immediately. ^ / a ul,u.,u, 

Theorem 2 The time average of I{x) converges in dis 
tribution: 



n-l 

■E 

k=0 



a.a.b 



(15) 



a = /3 - 1 

Proof. It is shown that the time average of the L]^ (m) 
functior0 converges to zero [l^, that is, for all / G 



where the random variable Ya^a,b has the following p.d.f, 

G'a.a.bi^) = 

(a — 6)sin7ra (a: — (a— a:) 



^ n— 1 

— f{T''x) -^0 as n ^ OO for almost all x. 



(21) 



(x-f))2a+2(2;-6)o(a-2:)o cos ira+(a-x)2o ('^ > ^) By thc thcorcm ([2T|) ■ for all e > there exist A^i such 

(16) that for n > A^i 



{b — a;)-^"+2(6— — a)*^ cos 7rQ + {a: — a)-^ 

where a ~ (3—1, that is, the random variables Ya^a,b 
obeys the generalized arcsine law. 

Proof. The time average of I{x) can be rewritten as 
lg/(TM=°^" + ^^"-^"^ (17) 



- n — 1 n — 1 

-V55,M(r'^x)-- V/(r'^x) 



fc=0 



n 



< - E |55,a,6(T'^2;) - /(T^x) 

Ti ^ ^ ' 



fc=0 



fc=0 



fc=0 



Using Theorem [I] we can write 

I ~ E < = Pr I (a 



<e/2, 
where 



6 < 



9S,a,b{.x) 



1 - G„ ( fr! 



(a > h) 

as n oo (18) 

(a < h). 

The derivative of (fTS)) . which is the p.d.f. of the time 
average of I{x) denoted by Ga,a,b, gives □ 

Definition 1 (i;\,c ,„ function v^fith finite mean) If the 

following conditions 



a for x G [0, 6) 
g{x) for X e [d, I - S] 
b forxe (1-5, 1]. 



(22) 



(23) 



Then by the continuity of g{x) and ((2T|) . for all e > 
there exist (5, and N2 such that for n > N2 



ji-i 



n-l 



i^5(T'^x)-1^5.-.,a,;,(r 



A:=0 



fe=0 



< e/2. 



(24) 



lim — < 00 

X dm 

and /or all e > 

f.l-£ 

Igjc?™ < 00 



(19) 



Therefore for all e > there exists (5* and N such 
that for n> N 



^ 71 — 1 



1 " " 

5(T^x) - - ^ /(T'^x) 



fc=0 



fe=0 



< e. 



(25) 



(20) By Theorem [51 the time average of g{x) converges in 



hold, then the function g is called the Lj^^{0, 1) function 
with finite mean with respect to m, denoted by L]^^ m('^' 1) 
function with finite mean. 

Examples of the L\^^{Q, 1) function w.r.t. the invari- 
ant measure of the modified Bernoulli map are I{x) and 
g{x) = x0 

^ Roughly speaking, we can say the L]„^^{0, 1) function 
with finite mean is considered as the L°°(0, 1) function. 



distribution to Yn 



□ 



Figures [U [2] and [3] demonstrate numerically that the 
p.d.f. for the time average of g{x) = x"^ obeys the theo- 
retical one even when the time n is finite {n = 10^). 

The following proof is not changed even if the condition 
of g{x) changes from the L\^^ m(Oi 1) function to the Z/°°(0, 1) 
function. 

* The L\(rn) function is the L^{m) function whose value 
is positive. 
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Fig. 1 The probability density function for the time average 
of g{x) =x^ [B^ 2.2). 



Fig. 3 The probability density function for the time average 
of g{x) ^x^ {B ^ 3.0). 




Fig. 2 The probability density function for the time average 
of g{x) = (B = 2.5). 



4 Application to the correlation functions 

We apply Theorem [3] to the correlation function, which 
is defined by the time average: 



(26) 



fc=0 



Corollary 1 For all n the correlation function of (7{x) 
converges to 1: 

1 

C{n) = lim — V a{xk)(j{xk+n) = 1, (27) 

N — ^Do /V < ^ 



k=0 



where a{x) = 1 (x £ [0,1/2)), -1 (x £ [1/2,1]). 

Proof. For all n the observation function gn (x) is defined 



gn{x) = a{x)aiT^x) 

' +1 a; e [0, a„] U [1 -a„,l] U A„ 
— 1 otherwise 



(28) 



where a„ = a„+i + 2^^^a^_^_l (a„ > and oq = 1/2) 
and An is the set which attains a{x)a-{T^'^x) = 1 and is 
subset of [a„, 1 — a„]: 



A„ ^ {x e [a„, 1 - a„] : a{x)a{T"x) = 1}. 



(29) 



Then .9ti(0) = gn{^) = 1 and gn{x) is continuous in 
[0, a„] U [1 — a„, 1] and for all e > 



gn{x)dm < I dm < oo. 



(30) 



In the case of finite measure, ergodic theory states that 
the correlation function defined by the time average equals 
to the correlation function defined by the ensemble av- 
erage, namely the average oi g{x)h{T"'x) with respect to 
the invariant measure. However, when the invariant mea- 
sure is not finite, the time-averaged correlation function 
is not equal to the ensemble-averaged correlation func- 
tionE In this section we demonstrate that the correlation 
function is intrinsically random in the modified Bernoulli 
map. 

^ The dependence of the decay of the ensemble-averaged 
correlation function on the initial ensemble is discussed in 



By Theorem [31 C (n) is convergence in distribution to 
Yct.i.i for all n. □ 

Corollary 2 For all n the correlation function of x is 
convergence in distribution to 5^0,0, i- 



C{n) = lim ^ XkXk+n ^ Yafl,!- 

N^OQ I\ ^ — ^ 



(31) 



Proof. For all n the observation function (/„ (x) is de- 
fined as 



gn{x) = x{T"x), 
and for all e > 

gn{x)dm < I dm < oo. 



(32) 
(33) 
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Then g„(0) — and gn{l) ~ 1, and gn{x) is continuous 
in [0, a„]U[l — a„, 1]. By Theoreni|3l C{n) is convergence 
in distribution to ^0,0,1 for all n. □ 

Figures [4] and [5] show that the correlation functions 
with fixed time difference n are intrinsically random and 
these distributions obey the generalized arcsine distribu- 
tion. 




Fig. 4 The probability density function for the correlation 
function defined by the time average of gio{x) = xxio {B — 
2.5). 



ft. 




0.4 0.6 

X 



Fig. 5 The probability density function for the correlation 
function defined by the time average of gio{x) = xxio {B — 
3.0). 



5 Stable law 

In this section we demonstrate the distribution for the 
time average of the Lj^^{0, 1) fimction with infinite mean, 
which satisfies the condition ([^0]) and 



lim ■ 

e— 



\g\dm 



00. 



(34) 



(35) 



In the case of the finite measure, the time average of the 
L[q^(0, 1) function with infinite mean converges to the 
stable distribution, as shown in Appendix A. 

In what follows, we study numerically the distribu- 
tion for the time average of the observation function 

x-" (x < i) 

(1-x-)"" ix>l), 

which is the Lj^^^O, 1) function with infinite mean and 
a > 0. As shown in Fig. 6, we find that the distribution 
for the scaled time average of Ll^^{Q, 1) function with in- 
finite mean also converges to the stable distribution with 
index 7. In numerical simulations we calculate the time 
average of g{x) for three different length of the simula- 
tion time (n = lO**, 10^ and 10^) and then determine the 
exponent 7 so as to make the distributions of the scaled 
time average invariant. Figure [7] shows the linear rela- 
tion between 7 and a clearly. Moreover, we find that the 
scaling exponent 7 obeys the non-trivial relation to the 
exponent a and B except for the case B < 2.5, i.e., 

7=;^ + l. (36) 

This relation is clearly shown in Fig. [5] except for the 
case B < 2.5.0 We summarize these results as following 
conjecture. 

Conjecture 1 Let g{x) he the L\^^{0, 1) function with 
infinite mean w.r.t. the invariant measure m. Further, 

x'^gix) = 0(1), .T ^ (37) 

il~xrg{x)^Oil), x^l (38) 

Then the scaled time average of g{x) converges to the 
stable distribution G^: 



n — 1 

where &„ oc n'^ and 7 



(39) 



B-l 



1. 



Remark 2 In the case of B ~ 2.0 the scaling sequence 
bn had better be n'^/logn rather than n'^ . 

Remark 3 The distribution for the scaled time average 
of g(x) = x^" also converges to the stable distribution 
when the invariant measure is finite. On the other hand, 
the distribution for the scaled time average of g{x) is dif- 
ferent from the stable one due to the generalized arcsine 
law for the occupation time of the interval [1/2, 1]. 

^ The reason that the relationship (|36p does not hold for 
the case B < 2.5 seems to be that the observation time is 
enoug in numerical simulations. 
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Fig. 8 Relationship between 7 and a and -B. 
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Fig. 6 The probability density function for the scaled time 
average of g{x) {B — 3.0 and a — 2.0). The fitting curve is a 
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Fig. 7 Scaling exponent 7 vs. a. 



6 Discussion 

In this paper we present the distributional limit theorems 
for the time average of the Lj^^ ^(0, 1) function with fi- 
nite mean and infinite mean using the modified Bernoulli 
map whose invariant measure is infinite. By applying this 
theorem to the correlation function, it is clearly shown 
that the correlation functions with fixed time difference n 
converge to the generalized arcsine distribution. G. Mar- 
golin and E. Barkai analyzed the distribution for the cor- 
relation function of a dichotomous random process and 
its convergence process [lOj . Our results correspond to 
the generalized result of their work in the way that the 
observation function can be extended to the Lj^^ ^(0, 1) 
function with finite mean^ However, the convergence 
process of the time average of the Lj^cTni^^ 1) ^^^^ finite 
mean is not studied. So we will study the convergence 



^ In [TD] the observation function is the characteristic func- 
tion I{x) with a = 1 and b = 0, which is the special case of 
our results. 



process for the correlation function in the forthcoming 
paper. 

From the viewpoint of physical observation, the dis- 
tributional limit theorems, namely, the Mittag-Leffler 
distribution and the generalized arcsine law and the sta- 
ble law, suggest that one can characterize the behaviour 
of non-stationary phenomena through the distribution 
of the time average. It is important to know what class 
the observation function belongs to, that is, whether the 
observation function is the Lj^^ ,n(0' 1) function with fi- 
nite mean or not. Because the observation function in 
physical systems is not always the L^{m) function. Ac- 
tually, we show that the correlation function is a typical 
example of the Lj^^ ^^{0,1) function with finite mean. 
This means that the correlation function, or the statis- 
tical quantities based on the time average in the inter- 
mittent phenomena, is intrinsically random. Universal 
distributions for the time average of various observation 
functions are shown in Table [1] 

In the concept of the "ergodicity" proposed by L. Boltz- 
mann, the ergodicity, i.e., the time average equals to the 
space average, guarantees the existence of the equilib- 
rium state in dynamical systems. In the non-equilibrium 
steady state, Sinai-Ruelle-Bowen measure is considered 
to describe the non-equilibrium steady state. However, 
there are no concepts of the "ergodicity" in the non-equilibrium 
state, that is, the measure of the non-equilibrium non- 
stationary state is not at all elucidated on the basis of 
the time average of the dynamical systems. El We hope 
that randomness of the time average in infinite mea- 
sure systems will give us the motive argument toward 
the ergodic problem in the non-equilibrium state. There 
is a possibility that Mittag-Leffler distribution or the 
generalized arcsine distribution or the stable distirbu- 
tion could become one of the measures characterizing 
the non-equilibrium state. Actually, these distribution 
universally appear in diffusion and its generalizations. 



From the aspect of the ensemble average, the approach to 
equilibrium of the o bservables is clearly shown using infinite 
measure systems in [21II22IP5] . 
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Table 1 Universal distributions of the time average of the observation function g{x). 



Invariant measure g(x) Distribution 



Finite {B < 2) 






Delta 


Finite {B < 2) 


L]^Jm) with infinite mean 


Stable 


Infinite {B > 2) 


L^{m) 
ri 




Mittag-Lejfler 


Infinite {B > 2) 


with finite mean 


Generalized arcsine 


Infinite {B > 2) 


ri 

loc.ra 


with infinite mean 


Stable 



Moreover, the generalized arcsine law has drawn much 
attention in disordered systems ^j,28j. 
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A The distribution for the time average of the 
non-i^(m) function in the case of the finite 
measure 

In the case of the finite measure (_B < 2), the invariant density 
can be written as 

p{x) = ^-^{x'-'' + {l~xf~^}. (40) 

Birkhoff 's ergodic theorem [29] tells us that the p.d.f. of the 
sequence {Tx, T^x, ■ ■ ■ , T'^x} obeys the invariant density as 
n —> oo. Let Xn be random variables with p.d.f. (|40p and 
g{x) = a;-" (a > 2 - B). The distribution of Y„ = ff(X„) is 
given by 

F(x) = Pr(Y <x) = Pr(X > x~i) 

= l-x'^. (41) 

The general central limit theorem [S^ says that the random 

variable (YiH 1- Y„)/n"/'^~^^ has the stable distribution. 

Therefore 

^ n—l 

^y^.9(r'°-)->Gf2-BWc asn^oo, (42) 

where G(2-B)/a is the stable distribution with index (2 — 
B)/a. 
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